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Abstract — Semiclassical models for multiple-user optical com- 
munication cannot assess the ultimate limits on reliable com- 
munication as permitted by the laws of physics. In all optical 
communications settings that have been analyzed within a quan- 
tum framework so far, the gaps between the quantum limit to the 
capacity and the Shannon limit for structured receivers become 
most significant in the low photon-number regime. Here, we 
present a quantum treatment of a multiple-transmitter multiple- 
receiver multi-spatial-mode free-space interference channel with 
diffraction-limited loss and a thermal background. We consider 
the performance of a laser-light (coherent state) encoding in 
conjunction with various detection strategies such as homodyne, 
heterodyne, and joint detection. Joint detection outperforms 
both homodyne and heterodyne detection whenever the channel 
exhibits "very strong" interference. We determine the capacity re- 
gion for homodyne or heterodyne detection when the channel has 
"strong" interference, and we conjecture the existence of a joint 
detection strategy that outperforms the former two strategies in 
this case. Finally, we determine the Han-Kobayashi achievable 
rate regions for both homodyne and heterodyne detection and 
compare them to a region achievable by a conjectured joint 
detection strategy. In these latter cases, we determine achievable 
rate regions if the receivers employ a recently discovered min- 
entropy quantum simultaneous decoder. 

I. Introduction 

The principal goals of information theory are to deter- 
mine the ultimate limits on reliable communication and to 
find ways of approaching these limits in practice. Point- 
to-point optical communication using laser-light modulation 
in conjunction with direct-detection and coherent-detection 
receivers has been studied in detail using the semiclassical 
theory of photodetection Q). This approach treats light as 
a classical electromagnetic field, and the fundamental noise 
encountered in photodetection is the shot noise associated with 
the discreteness of the electron charge. 

These semiclassical treatments for systems that exploit 
classical-light modulation and conventional receivers (direct, 
homodyne, or heterodyne) have had some success, but we 
should recall that electromagnetic waves are quantized, and 
the correct assessment of systems that use non-classical light 
sources and/or general optical measurements requires a full 
quantum-mechanical framework |2|. Consider several recent 
theoretical studies on the point-to-point |3], [4|, broadcast 
and multiple-access [6| bosonic channels with linear loss and 
a thermal background. These studies have shown that achiev- 
able communication rates surpass what can be obtained with 
conventional receivers. Prior work has established that Holevo 



information rates are achievable for information transmission 
on general quantum point-to-point J7), (8), broadcast ||9ll 
and multiple-access [10| channels. In each case, the Holevo 
information rates are an upper bound to the Shannon rates 
computed for any specific transmitter-modulation receiver- 
measurement pair. For the general quantum channel, attaining 
Holevo information rates may require collective measurements 
(a joint detection) across the channel outputs. 

The next level of complexity beyond the point-to-point, 
broadcast, or multiple-access channels is arguably captured 
by the interference channel ifTTIl . fl2l . fl3l . This channel, in 
general, can have M senders and M receivers (M > 2), where 
each sender would like to communicate only with a partner 
receiver^] but most research focuses on the special case of two 
senders and two receivers. The Gaussian interference channel 
has been analyzed in depth in the classical information-theory 
literature, and this model readily applies for an optical interfer- 
ence channel with coherent-state inputs and coherent detection. 
Calculating the capacity of this channel in the general case has 
been an open problem for some time, but several researchers 
have found it in the special cases of "very strong" [ 1 1 1 and 
"strong" interference lfl2ll . 1131 . Also, Han and Kobayashi 
determined the best inner bound on the channel's capacity, 
by having each receiver partially decode the message of the 
other sender along with a full decoding of the message of the 
partner sender fil l. However, all of these strategies assume 
that the underlying channel is classical, and we would expect 
a quantum strategy with a power-constrained encoding and 
collective measurement at the receivers to outperform such 
strategies. 

In this paper, we consider a pure-loss thermal-noise bosonic 
interference channel, particularly in the context of free- 
space (wireless) terrestrial optical communications. We as- 
sume a coherent-state encoding throughout this paper. We find 
achievable rate regions for the two-sender two-receiver chan- 
nel with coherent-detection receivers, for the "very strong" 
and "strong" interference regimes. We find the rate region 
achievable with a joint-detection receiver (JDR) under "very 
strong" interference. We also determine a "min-entropy" JDR- 
achievable rate region for the "strong" interference regime 

Note that this interference channel setting has been extensively analyzed 
in network information theory, for scaling behavior of the total source-to- 
destination capacity and communication latency as M grows, for randomly 
distributed M source-destination node pairs in a given network area. 



by exploiting a recent result from Ref. |[T4"1 . and we find a 
conjectured JDR-achievable region were a conjecture from 
Ref. [14] regarding quantum simultaneous decoding true (c.f, 
page 4-15 of Ref. [15] for a classical simultaneous decoder). 
Next, we evaluate the Han-Kobayashi rate region for ho- 
modyne and heterodyne detection in the general case, and 
we show an achievable rate region using a "min-entropy" 
quantum simultaneous decoder from Ref. Ifl4l . Finally, we 
conjecture a Han-Kobayashi-like achievable rate region using 
Conjecture 2 from Ref. [14]. Our results here differ from those 
in Ref. lfT4l — there some of us considered the general quantum 
interference channel whereas here we consider specifically a 
bosonic interference channel. 

II. A Free-Space Optical Interference Channel 

Consider a range-i line-of-sight free-space optical channel 
with hard circular transmit and receive apertures of areas 
A t and A r respectively. Assume A-center-wavelength quasi- 
monochromatic transmission. In the near-field propagation 
regime (Fresnel number product, Df = A t A r /(XL) 2 ^> 1), a 
normal-mode decomposition of the free-space optical channel 
yields M ~ 2Df orthogonal spatio-polarization transmitter- 
to-receiver modes (D t spatial modes, each of two orthogonal 
polarizations) with near-unity transmitter-to-receiver power 
transmissivities {rj m w 1). In the far-field propagation regime 
(Df <C 1), only two orthogonal spatial modes (one of each 
orthogonal polarization) have appreciable power transmissivity 
(77 Df for each mode). 

Sender m modulates her information on the m th transmitter- 
pupil spatial mode, and Receiver m separates and demodu- 
lates information from the corresponding receiver-pupil spatial 
mode. With perfect spatial-mode control at the transmitter and 
perfect mode separation at the receiver, the orthogonal spatial 
modes can be thought of as independent parallel channels with 
no cross talk. However, imperfect (slightly non-orthogonal) 
mode generation or imperfect mode separation can result in 
cross talk (interference) between the M channels. 

We take our interference channel model as a passive linear 
mixing of the input modes along with the possibility of a 
thermal environment adding zero-mean, isotropic Gaussian 
noise. Although the results here apply to cyclic interference 
channels with M senders and M receivers, for simplicity, we 
limit ourselves to the M = 2 case, in which case the channel 
model reduces to 

h = y/rjiiai + y/rj^xa-z + VmPi, (1) 
h = y/^&i - y/rj^a 2 + y%i>2, (2) 

where 7711, 7712, 7721,7722,771,772 € y/VuJfri = ^2 17722, 
f\x = l-r\xx - 7721, and fj 2 = 1 - J7 12 - 7722- The following 
conditions ensure that the network is passive: 

»7i:L+J7i2 < 1) rin+mi < 1, 7722+7721 < 1, 7722+7712 < 1- 

We constrain the mean photon number of the transmitters Si 
and a,2 to be Ns ± and Ns 2 photons per mode, respectively, the 
environment modes v\ and O2 are in statistically independent 



zero-mean thermal states ||2| with respective mean photon 
numbers Nb 1 and Nb 2 per mode. 

Note that for a coherent-state encoding and coherent- 
detection at both receivers, the above model is a special case 
of a complex Gaussian interference channel, and we can study 
its capacity regions in various settings by applying the results 
from Refs. ATI . HV2\ , ifTJl . If the senders prepare their inputs 
in coherent states \a>i) and (0:2), with a\,a2 G K, and 
both receivers perform real-quadrature homodyne detection 
on their respective modes, the result is a classical Gaussian 
interference channel j2), where Receivers 1 and 2 obtain 
respective conditional Gaussian random variables Yj and Y 2 
distributed as 

Y\ ~ JV(y77iTai + % /772Ta 2 , (2fh.N Bl + 1) /4) , 
Y 2 ~N (0712^2 + y/f^a u {2f) 2 N B2 + 1) /4) , 

where the term in the noise variances arises physically 

from the zero-point fluctuations of the vacuum. Suppose that 
the senders again encode their signals as coherent states \a\) 
and \a 2 ), but this time with ot\, ct2 € C, and that the receivers 
both perform heterodyne detection. This results in a classical 
complex Gaussian interference channel [2], where Receivers 1 
and 2 detect respective conditional complex Gaussian random 
variables Z x and Z 2 , whose real parts are distributed as 

Re {Z m } ~ TV (/j m , {fj m N Bm + l)/2) , 

where m G {1,2}, /Ji = ^/rj^Re {ai} + v /?72iRe {a 2 }, 
/i 2 = Y / 77i2 Re {cti} + -y/7722 Re {a 2 }, and the imaginary parts 
of Z\ and Z 2 are distributed with the same variance as their 
real parts, and their respective means are ^7711 Im {a\} + 
■^7721 Im {«2} and ^7712 Im {a±} + ^7722 Im {02}. The factor 
of 1/2 in the noise variances is due to the attempt to measure 
both quadratures of the field simultaneously [2|. 

III. Very Strong Interference 

Carleial determined the capacity region of a classical Gaus- 
sian interference channel in the case of "very strong" inter- 
ference ifTTl . Suppose that X\ and X 2 are the input random 
variables for Senders 1 and 2 and that B\ and B2 represent 
the outputs for Receivers 1 and 2, respectively. Then the most 
general way to state the condition for very strong interference 
is that the following information inequalities should hold for 
all input distributions pxi(+i) and px 2 ( x 2) 021 : 

I{X 1 ;B 1 \X 2 )<I{X 1 -B 2 ), (3) 
I(X 2 ;B 2 \X 1 )<I{X 2 ;B 1 ). (4) 

Carleial proved the surprising result [11] that the capacity 
region of the interference channel under this setting is the 
convex closure of positive rate pairs (Ri,R 2 ) such that 

iii</(Jfi;Bi|Jf 2 ), R2<I{X2;B 2 \X 1 ), (5) 

for some input distributions px ± (#1) and px 2 (x 2 ). The coding 
strategy is simply for each receiver to decode first what the 
other sender is transmitting, remove this signal, and then 
decode the message intended for him. 



The conditions in (3]4i translate to the following ones for 
the case of coherent-state encoding and coherent detection: 
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where i = 1 for homodyne detection and i = for heterodyne 
detection. 

We can also consider the case when the senders employ 
coherent-state encodings and the receivers employ a joint- 
detection strategy on all of their respective channel outputs. 
The conditions in ( 3]|4 1 readily translate to this quantum setting 
where we now consider B\ and B 2 to be quantum systems, 
and the information quantities in ([3]|4]i and |5]l now become 
Holevo informations fl4l . Theorem 6 in Ref. 01411 provides 
a simple proof that Carleial's result applies in the quantum 
domain to these Holevo informations. So, the conditions in 
djpl) when restricting to coherent-state encodings translate to 



9 {V22Ns 2 + fj 2 N B2 ) - g {mN B2 ) (8) 

< 9 (V2iN S2 + VuN Sl + ViN Bl ) - 9 (VnN Sl + ViN Bl ) , 

9 {miN Sl + mN Bl ) - 9 imN Bl ) (9) 

< g (r)i2N Sl + V22N S2 + f) 2 N B 2 ) - 9 (V22N S2 + m N B 2 ) ■ 

where g (N) is the entropy of a thermal state with mean photon 
number N: 

g (N) = (N + 1) In (N + 1) - N In (N) . 

An achievable rate region is then 

Ri < g{miN Sl +fjiNB 1 )-g(fjiN Bl ), 
R2 < g {V22Ns 2 + mN B2 ) - g (fj 2 N B2 ) . 

These rates are achievable using a coherent-state encoding, 
but not necessarily optimal (though they would be optimal if 
the minimum-output entropy conjecture from Refs. J5], |[T6ll 
were true). Nevertheless, these rates always beat the rates from 
homodyne and heterodyne detection, and Figure [T] displays 
examples of the capacity (and achievable rate) regions in the 
low- and high-power regimes. 

IV. Strong Interference 

Sato 1 12 1 and Han-Kobayashi lfl3l independently deter- 
mined the capacity of a classical Gaussian interference channel 
under "strong" interference. A channel has "strong" inter- 
ference if the following information inequalities hold for all 

PXii^i) and px 2 (x 2 ) US): 



I(X 1 ;B 1 \X 2 )<I(X 1 ;B 2 \X 2 ) 
I{X 2 ;B 2 \X 1 )<I(X 2 ;B 1 \X 1 ) 



(10) 
(11) 
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Fig. 1 . Capacity regions for coherent-state encodings and coherent detection, 
and achievable rate regions for coherent-state encodings and joint-detection 
receivers — both with 7711 = -(722 = 1/16 and 7712 = ?72i = 1/2 ("very 
strong" interference for coherent detection and such that \&\9) hold). The 
LHS displays these regions in a low-power regime with Ns x = N$ 2 = 1 
and Ng 1 = Ng 2 = 1, and the RHS displays these regions in a high- 
power regime where Ng 1 = Ng 2 = 100. Homodyne detection outperforms 
heterodyne detection in the low-power regime because it has a reduced 
detection noise, while heterodyne detection outperforms homodyne detection 
in the high-power regime because its has an increased bandwidth. 



The capacity region of the classical interference channel under 
this setting is the convex closure of positive rate pairs (i2j , R 2 ) 
such that E), US, 03 : 



Rx <I{Xx;B i \X 2 ) ) R2<I(X 2 ;B 2 \X 1 ), 
R x +R 2 < ram{I {XxX^Bi) ,1 (X X X 2 ;B 2 )} . 



(12) 
(13) 



The conditions in ( flOp"!) translate to the following ones for 
coherent-state encoding and coherent detection: 
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and the capacity region has the two inequalities in (6]|7i and 
an additional bound on the sum rate: 
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where again i = 1 for homodyne detection and i 
heterodyne detection. 

The situation for a joint-detection strategy over all of the 
channel outputs becomes more complicated for the case of 
"strong" interference, because we require a quantum simulta- 
neous decoder [14] in order to achieve the information rates in 
( 12p3 i with B\ and B 2 becoming quantum systems. Such a 
simultaneous decoder is analogous to a classical simultaneous 
decoder (e.g., see page 4-15 of Ref. [15]), but we have not 
yet been able to prove the existence of it in the quantum case 
(see Conjecture 2 of Ref. lfT4l ). Yet, we do have an achievable 
simultaneous decoding strategy expressed in terms of min- 
entropies (see Theorem 4 of Ref. 1 14 1), where the min-entropy 
of a probability distribution is the negative logarithm of the 
probability of its mode ifTTI . and, as a simple extension of 
this idea, the min-entropy of a density operator is the negative 
logarithm of its maximum eigenvalue. For a thermal state with 
average photon number N B , its min-entropy is In (N B + 1), 
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Fig. 2. The parameters of the free-space channel are Ng 1 = Ng 2 = 1, 
771 1 = »722 = 0.3, T)2i = r?i2 = 0.6, Ng 1 = Ng 2 = 2 for the low-power 
regime, and Ns t = Ns 2 = 100 for the high-power regime. The LHS figure 
depicts the "strong" interference capacity regions in the low-power regime 
for homodyne and heterodyne detection, the achievable rate region with the 
conjectured joint detection, and the convex hull over regions arising from 
different min-entropy decoding strategies. Homodyne detection outperforms 
heterodyne detection in this regime, and the min-entropy decoders perform the 
worst. The RHS figure displays the same regions in the high-power regime. 
Heterodyne detection outperforms homodyne detection, and interestingly, the 
min-entropy decoders outperform heterodyne detection. 



and this result allows us to determine an achievable rate 
region with a simultaneous decoding strategy (similar to the 
simultaneous decoding inner bound on page 6-7 of Ref. 1 15 1): 

Ri < In (miN Sl + mN Bl +l)-g (fjiN Bl ) , 
R 2 < In (?722iVs 2 + f] 2 N B2 +l)-g (fj 2 N B2 ) , 

Ri + R2 < 

g (viiNsi + mi N s 2 + viNbi) - g {fjiN Bl ) , 
g (r)22Ns 2 + m2N Sl + mN B2 ) - g (fj 2 N B2 ) 

This is one particular variation of an achievable strategy in 
which we have the bounds on the individual rates expressed 
in terms of min-entropies and the bound on the sum rate 
expressed with von Neumann entropies, though note that there 
are other variations we could consider in light of Theorem 4 of 
Ref. [14|. We can then take the convex hull of the achievable 
rate regions for these different strategies to get an achievable 
rate region for a min-entropy quantum simultaneous decoder 
(the RHS of Figure [2] displays an interesting example of such a 
"min-entropy" region). If Conjecture 2 of Ref. fPfll regarding 
the existence of a quantum simultaneous decoder were true, 
then the rate region in (|12||13| l would be achievable under 
under the conditions of < | 1 OH 1 1 j > (with Holevo information rates 
replacing Shannon rates). Figure|2]displays the different capac- 
ity and achievable rate regions when a free-space interference 
channel exhibits "strong" interference. 

V. Han-Kobayashi Rate Regions 

The Han-Kobayashi region is the best known achievable rate 
region for the classical interference channel fl3l . The coding 
strategy to achieve this region is for each receiver to decode 
partially the other sender's message while fully decoding the 
partner sender's message. With this strategy, the four parties 
can choose to take advantage of channel interference while 
achieving the task of paired sender-receiver communication. 



A compact description of the Han-Kobayashi region comes 
from its reduction with a Fourier-Motzkin elimination al- 
gorithm [18). It is the convex closure of all positive rate 
pairs (Ri,R 2 ) satisfying the following inequalities and the 
inequalities obtained from the ones below by swapping the 
indices 1 and 2: 

iZi</(EWi;Bi|W2), (14) 

Ri < / [U\] Bi\W\W 2 ) + I (Wr, B 2 \U 2 W 2 ) , (15) 

R1 + R2 <I{U 1 ;B 1 \W 1 W 2 )+I{U 2 W 2 W 1 ;R 2 ), (16) 

R1+R2 <I{UiW 2l B l \W 1 )+I{U 2 W 1 -B 2 \W 2 ) 1 (17) 
2Ri+R 2 <I(Ui;B 1 \WiW 2 ) + I(U 1 W 1 W 2 ;B 1 ) 

+ I{U 2 W 1 ;B 2 \W 2 ). (18) 

In the above, U m is the "personal" random variable of 
Sender m, and W m is her "common" random variable. 

The Han-Kobayashi coding strategy readily translates into 
a strategy for coherent-state encoding and coherent detection. 
Sender m shares the total photon number Ng m between her 
personal message and her common message. Let A m be the 
fraction of signal power that Sender m devotes to her personal 
message, and let A m denote the other fraction of signal 
power that Sender m devotes to her common message. The 
inequalities above become the following ones for the case of 
coherent-state encoding along with coherent detection: 
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where 



7 (x) = In (1 + x) /2\ N m = [2% n N Bm + l] /A\ 

i = 1 for homodyne detection, i = for heterodyne detection, 
and m <E {1, 2}. 

We also conjecture a Han-Kobayashi achievable rate re- 
gion if the senders employ coherent-state encodings and the 
receivers exploit joint-detection receivers (this again follows 
from Conjecture 2 of Ref. [14| regarding the existence of 
a quantum simultaneous decoder). The inequalities for the 
region are similar to those in ( T4p8 1 and the additional 



Han-Kobayashi rate region for different detectors 




Ri 

Fig. 3. The free-space channel no longer exhibits "very strong" or "strong" 
interference for coherent detection because the parameters are N,s 1 = Ns 2 = 
100, N Bl = Nb 2 = 1. Wu = V22 = 0.8, and r/2i = r?i 2 = 0.1. The figure 
depicts the achievable rate regions by employing a Han-Kobayashi coding 
strategy for homodyne and heterodyne detection. It also depicts what we can 
achieve by employing 49 variations of the min-entropy decoding strategy from 
Ref. 1 14] and taking the convex hull of these achievable rate regions. Finally, 
it displays the conjectured region if a joint detection decoder were to exist 
(see Conjecture 2 of Ref. 1141 ). All of these regions are with respect to a 
10%-personal, 90%-common Han-Kobayashi power split. 

"swapped" inequalities, with the exception that Holevo infor- 
mations replace mutual informations. 

We can obtain an achievable rate region by exploiting the 
quantum simultaneous decoder from Theorem 4 of Ref. [14] 
that gives rates which are a difference of a min-entropy 
and a von Neumann entropy. The region's characterization 
is in terms of the Han-Kobayashi (HK) characterization with 
14 inequalities lfl3ll , corresponding to two different multiple 
access channels (MACs) induced to each receiver by the HK 
coding strategy (seven inequalities for each MAC). There 
are 49 variations of these min-entropy decoders — one of the 
information rates in the seven inequalities for each MAC is von 
Neumann and the other six are min-entropy rates. Then taking 
the convex hull of these 49 different achievable rate regions 
gives an achievable rate region for a min-entropy decoding 
strategy. Figure [3] plots the regions achievable with coherent 
detection, the min-entropy decoder, and the conjectured joint 
detector for a particular HK power split. 

VI. Conclusion 
The semiclassical models for free-space optical commu- 
nication are not sufficient to understand the ultimate limits 
on reliable communication rates, for both point-to-point and 
multiple-sender-receiver channels. We presented a quantum- 
mechanical model for the free-space optical interference 
channel and determined achievable rate regions using both 
structured and unstructured receivers. Interestingly, the min- 
entropy decoder from Ref. Ifl4l can achieve rates that are 
unachievable by both homodyne and heterodyne detection 



when the channel exhibits "strong interference." Finally, we 
determined the Han-Kobayashi inner bound for homodyne and 
heterodyne detection, and we conjectured a rate region of this 
form if a quantum simultaneous decoder were to exist. 

Several open problems remain for this line of inquiry. Per- 
haps the biggest open question is to prove Conjecture 2 from 
Ref. [14 1 concerning the existence of a quantum simultaneous 
decoder for a general quantum interference channel. Also, we 
do not know if a coherent-state encoding is in fact optimal for 
the free-space interference channel — it might be that squeezed 
state transmitters could achieve higher communication rates as 
in Ref. f6|. One could also evaluate the ergodic and outage 
capacity regions based on the statistics of rjij, which could be 
derived from the spatial coherence functions of the stochastic 
mode patterns under atmospheric turbulence. 
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